Flexibility in the Periplasmic Domain of BamA Is Important for Function.

The ß-barrel assembly machine (BAM) mediates the biogenesis of outer membrane proteins (OMPs) in Gram-negative bacteria. BamA, the central BAM subunit composed of a transmembrane ß-barrel domain linked to five polypeptide transport-associated (POTRA) periplasmic domains, is thought to bind nascent OMPs and undergo conformational cycling to catalyze OMP folding and insertion. One model is that conformational flexibility between POTRA domains is part of this conformational cycling. Nuclear magnetic resonance (NMR) spectroscopy was used here to study the flexibility of the POTRA domains 1-5 in solution. NMR relaxation studies defined effective rotational correlational times and together with residual dipolar coupling data showed that POTRA1-2 is flexibly linked to POTRA3-5. Mutants of BamA that restrict flexibility between POTRA2 and POTRA3 by disulfide crosslinking displayed impaired function in vivo. Together these data strongly support a model in which conformational cycling of hinge motions between POTRA2 and POTRA3 in BamA is required for biological function.

Structure-Based Assignment of Ile, Leu, and Val Methyl Groups in the Active and Inactive Forms of the Mitogen-Activated Protein Kinase Extracellular Signal-Regulated Kinase 2.

Resonance assignments are the first step in most NMR studies of protein structure, function, and dynamics. Standard protein assignment methods employ through-bond backbone experiments on uniformly (13)C/(15)N-labeled proteins. For larger proteins, this through-bond assignment procedure often breaks down due to rapid relaxation and spectral overlap. The challenges involved in studies of larger proteins led to efficient methods for (13)C labeling of side chain methyl groups, which have favorable relaxation properties and high signal-to-noise. These methyls are often still assigned by linking them to the previously assigned backbone, thus limiting the applications for larger proteins. Here, a structure-based procedure is described for assignment of (13)C(1)H3-labeled methyls by comparing distance information obtained from three-dimensional methyl-methyl nuclear Overhauser effect (NOE) spectroscopy with the X-ray structure. The Ile, Leu, or Val (ILV) methyl type is determined by through-bond experiments, and the methyl-methyl NOE data are analyzed in combination with the known structure. A hierarchical approach was employed that maps the largest observed "NOE-methyl cluster" onto the structure. The combination of identification of ILV methyl type with mapping of the NOE-methyl clusters greatly simplifies the assignment process. This method was applied to the inactive and active forms of the 42-kDa ILV (13)C(1)H3-methyl labeled extracellular signal-regulated kinase 2 (ERK2), leading to assignment of 60% of the methyls, including 90% of Ile residues. A series of ILV to Ala mutants were analyzed, which helped confirm the assignments. These assignments were used to probe the local and long-range effects of ligand binding to inactive and active ERK2.

Dynamics of protein kinases: insights from nuclear magnetic resonance.

Protein kinases are ubiquitous enzymes with critical roles in cellular processes and pathology. As a result, researchers have studied their activity and regulatory mechanisms extensively. Thousands of X-ray structures give snapshots of the architectures of protein kinases in various states of activation and ligand binding. However, the extent of and manner by which protein motions and conformational dynamics underlie the function and regulation of these important enzymes is not well understood. Nuclear magnetic resonance (NMR) methods provide complementary information about protein conformation and dynamics in solution. However, until recently, the large size of these enzymes prevented researchers from using these methods with kinases. Developments in transverse relaxation-optimized spectroscopy (TROSY)-based techniques and more efficient isotope labeling strategies are now allowing researchers to carry out NMR studies on full-length protein kinases. In this Account, we describe recent insights into the role of dynamics in protein kinase regulation and catalysis that have been gained from NMR measurements of chemical shift changes and line broadening, residual dipolar couplings, and relaxation. These findings show strong associations between protein motion and events that control kinase activity. Dynamic and conformational changes occurring at ligand binding sites and other regulatory domains of these proteins propagate to conserved kinase core regions that mediate catalytic function. NMR measurements of slow time scale (microsecond to millisecond) motions also reveal that kinases carry out global exchange processes that synchronize multiple residues and allosteric interconversion between conformational states. Activating covalent modifications or ligand binding to form the Michaelis complex can induce these global processes. Inhibitors can also exploit the exchange properties of kinases by using conformational selection to form dynamically quenched states. These investigations have revealed that kinases are highly dynamic enzymes, whose regulation by interdomain interactions, ligand binding, and covalent modifications involve changes in motion and conformational equilibrium in a manner that can be correlated with function. Thus, NMR provides a unique window into the role of protein dynamics in kinase regulation and catalysis with important implications for drug design.

Slow inhibition and conformation selective properties of extracellular signal-regulated kinase 1 and 2 inhibitors.

The mitogen-activated protein (MAP) kinase pathway is a target for anticancer therapy, validated using inhibitors of B-Raf and MAP kinase kinase (MKK) 1 and 2. Clinical outcomes show a high frequency of acquired resistance in patient tumors, involving upregulation of activity of the MAP kinase, extracellular signal-regulated kinase (ERK) 1 and 2. Thus, inhibitors for ERK1/2 are potentially important for targeted therapeutics against cancer. The structures and potencies of different ERK inhibitors have been published, but their kinetic mechanisms have not been characterized. Here we perform enzyme kinetic studies on six representative ERK inhibitors, with potencies varying from 100 pM to 20 µM. Compounds with significant biological activity (IC50 < 100 nM) that inhibit in the subnanomolar range (Vertex-11e and SCH772984) display slow-onset inhibition and represent the first inhibitors of ERK2 known to demonstrate slow dissociation rate constants (values of 0.2 and 1.1 h(-1), respectively). Furthermore, we demonstrate using kinetic competition assays that Vertex-11e binds with differing affinities to ERK2 in its inactive, unphosphorylated and active, phosphorylated forms. Finally, two-dimensional heteronuclear multiple-quantum correlation nuclear magnetic resonance experiments reveal that distinct conformational states are formed in complexes of Vertex-11e with inactive and active ERK2. Importantly, two conformers interconvert in equilibrium in the active ERK2 apoenzyme, but Vertex-11e strongly shifts the equilibrium completely to one conformer. Thus, a high-affinity, slow dissociation inhibitor stabilizes different enzyme conformations depending on the activity state of ERK2 and reveals properties of conformational selection toward the active kinase.

Optimized reverse micelle surfactant system for high-resolution NMR spectroscopy of encapsulated proteins and nucleic acids dissolved in low viscosity fluids.

An optimized reverse micelle surfactant system has been developed for solution nuclear magnetic resonance studies of encapsulated proteins and nucleic acids dissolved in low viscosity fluids. Comprising the nonionic 1-decanoyl-rac-glycerol and the zwitterionic lauryldimethylamine-N-oxide (10MAG/LDAO), this mixture is shown to efficiently encapsulate a diverse set of proteins and nucleic acids. Chemical shift analyses of these systems show that high structural fidelity is achieved upon encapsulation. The 10MAG/LDAO surfactant system reduces the molecular reorientation time for encapsulated macromolecules larger than ~20 kDa leading to improved overall NMR performance. The 10MAG/LDAO system can also be used for solution NMR studies of lipid-modified proteins. New and efficient strategies for optimization of encapsulation conditions are described. 10MAG/LDAO performs well in both the low viscosity pentane and ultralow viscosity liquid ethane and therefore will serve as a general surfactant system for initiating solution NMR studies of proteins and nucleic acids.

Phosphorylation releases constraints to domain motion in ERK2.

Protein motions control enzyme catalysis through mechanisms that are incompletely understood. Here NMR (13)C relaxation dispersion experiments were used to monitor changes in side-chain motions that occur in response to activation by phosphorylation of the MAP kinase ERK2. NMR data for the methyl side chains on Ile, Leu, and Val residues showed changes in conformational exchange dynamics in the microsecond-to-millisecond time regime between the different activity states of ERK2. In inactive, unphosphorylated ERK2, localized conformational exchange was observed among methyl side chains, with little evidence for coupling between residues. Upon dual phosphorylation by MAP kinase kinase 1, the dynamics of assigned methyls in ERK2 were altered throughout the conserved kinase core, including many residues in the catalytic pocket. The majority of residues in active ERK2 fit to a single conformational exchange process, with kex ≈ 300 s(-1) (kAB ≈ 240 s(-1)/kBA ≈ 60 s(-1)) and pA/pB ≈ 20%/80%, suggesting global domain motions involving interconversion between two states. A mutant of ERK2, engineered to enhance conformational mobility at the hinge region linking the N- and C-terminal domains, also induced two-state conformational exchange throughout the kinase core, with exchange properties of kex ≈ 500 s(-1) (kAB ≈ 15 s(-1)/kBA ≈ 485 s(-1)) and pA/pB ≈ 97%/3%. Thus, phosphorylation and activation of ERK2 lead to a dramatic shift in conformational exchange dynamics, likely through release of constraints at the hinge.

Efficient detection of hydrogen bonds in dynamic regions of RNA by sensitivity-optimized NMR pulse sequences.

Improved Sensitivity: Efficient NMR experiments are presented for determining the secondary structure in large and dynamic RNAs using J-couplings across hydrogen bonds. The experiments provide up to eight-fold improved sensitivity and thus enable detection of base pairs in dynamic regions even in large RNAs.

Structure of the BamC two-domain protein obtained by Rosetta with a limited NMR data set.

The CS-RDC-NOE Rosetta program was used to generate the solution structure of a 27-kDa fragment of the Escherichia coli BamC protein from a limited set of NMR data. The BamC protein is a component of the essential five-protein ß-barrel assembly machine in E. coli. The first 100 residues in BamC were disordered in solution. The Rosetta calculations showed that BamC101₋344 forms two well-defined domains connected by an ~18-residue linker, where the relative orientation of the domains was not defined. Both domains adopt a helix-grip fold previously observed in the Bet v 1 superfamily. ¹5N relaxation data indicated a high degree of conformational flexibility for the linker connecting the N-terminal domain and the C-terminal domain in BamC. The results here show that CS-RDC-NOE Rosetta is robust and has a high tolerance for misassigned nuclear Overhauser effect restraints, greatly simplifying NMR structure determinations.

Structure and flexibility of the complete periplasmic domain of BamA: the protein insertion machine of the outer membrane.

Folding and insertion of ß-barrel outer membrane proteins (OMPs) is essential for Gram-negative bacteria. This process is mediated by the multiprotein complex BAM, composed of the essential ß-barrel OMP BamA and four lipoproteins (BamBCDE). The periplasmic domain of BamA is key for its function and contains five "polypeptide transport-associated" (POTRA) repeats. Here, we report the crystal structure of the POTRA4-5 tandem, containing the essential for BAM complex formation and cell viability POTRA5. The domain orientation observed in the crystal is validated by solution NMR and SAXS. Using previously determined structures of BamA POTRA1-4, we present a spliced model of the entire BamA periplasmic domain validated by SAXS. Solution scattering shows that conformational flexibility between POTRA2 and 3 gives rise to compact and extended conformations. The length of BamA in its extended conformation suggests that the protein may bridge the inner and outer membranes across the periplasmic space.

Ultrafast dynamics show that the theophylline and 3-methylxanthine aptamers employ a conformational capture mechanism for binding their ligands.

RNAs often exhibit a high degree of conformational dynamics and heterogeneity, leading to a rugged energy landscape. However, the roles of conformational heterogeneity and rapid dynamics in molecular recognition or RNA function have not been extensively elucidated. Ultrafast time-resolved fluorescence spectroscopic experiments were used here to probe picosecond dynamics of the theophylline-binding RNA aptamer. These studies showed that multiple conformations are populated in the free RNA, indicating that this aptamer employs a conformational capture mechanism for ligand binding. The base on residue 27 in an internal loop exists in at least three conformational states in the free RNA, including binding competent and incompetent states that have distinct fluorescence decay signatures indicating different base stacking interactions. Picosecond dynamics were also detected by anisotropy experiments, where these motions indicate additional dynamics for base 27. The picosecond data show that theophylline binding shifts the equilibrium for conformations of base 27 from primarily stacked in the free RNA to mostly unstacked in the RNA-theophylline complex, as observed in the previous NMR structure. In contrast, base 10 in a second internal loop is mostly preorganized in the free RNA, consistent with it being stacked between G11 and G25, as is observed in the bound state. Picosecond dynamics were also measured on a modified aptamer that binds with higher affinity to 3-methylxanthine than theophylline. The modified aptamer shows less heterogeneity in the aptamer-3-methylxanthine complex than what is observed in the theophylline aptamer-theophylline complex.

Longitudinal-relaxation-enhanced NMR experiments for the study of nucleic acids in solution.

Atomic-resolution information on the structure and dynamics of nucleic acids is essential for a better understanding of the mechanistic basis of many cellular processes. NMR spectroscopy is a powerful method for studying the structure and dynamics of nucleic acids; however, solution NMR studies are currently limited to relatively small nucleic acids at high concentrations. Thus, technological and methodological improvements that increase the experimental sensitivity and spectral resolution of NMR spectroscopy are required for studies of larger nucleic acids or protein-nucleic acid complexes. Here we introduce a series of imino-proton-detected NMR experiments that yield an over 2-fold increase in sensitivity compared to conventional pulse schemes. These methods can be applied to the detection of base pair interactions, RNA-ligand titration experiments, measurement of residual dipolar (15)N-(1)H couplings, and direct measurements of conformational transitions. These NMR experiments employ longitudinal spin relaxation enhancement techniques that have proven useful in protein NMR spectroscopy. The performance of these new experiments is demonstrated for a 10 kDa TAR-TAR*(GA) RNA kissing complex and a 26 kDa tRNA.

Chemical shift anisotropy of imino 15N nuclei in Watson-Crick base pairs from magic angle spinning liquid crystal NMR and nuclear spin relaxation.

Knowledge of (15)N chemical shift anisotropy is prerequisite both for quantitative interpretation of nuclear spin relaxation rates in terms of local dynamics and for the use of residual chemical shift anisotropy (RCSA) as a constraint in structure determination. Accurate measurement of the very small RCSA from the difference in (15)N chemical shift under isotropic and weakly aligning liquid crystalline conditions is very sensitive to minute differences in sample conditions, such as pH or ionic strength. For this reason, chemical shifts were measured for the same solution, under static liquid crystalline alignment, and under magic angle spinning conditions where alignment relative to the magnetic field is removed. Measurements were made for 14 well-resolved G-N(1) and 6 U-N(3) (15)N nuclei in a sample of tRNA(Val). Fitting these RCSA data together with (15)N-(1)H dipole-CSA cross-correlated relaxation measurements to the recently refined structural model of tRNA(Val) yields the magnitude, asymmetry, and orientation of the (15)N CSA tensors.

NMR chemical exchange as a probe for ligand-binding kinetics in a theophylline-binding RNA aptamer.

The apparent on and off rate constants for binding of theophylline to its RNA aptamer in the absence of Mg(2+) were determined here by 2D (1)H-(1)H ZZ-exchange NMR spectroscopy. Analysis of the buildup rate of the exchange cross peaks for several base-paired imino protons in the RNA yielded an apparent k(on) of 600 M(-1) s(-1). This small apparent k(on) results because the free RNA exist as a dynamic equilibrium of inactive states rapidly interconverting with a low population of active species. The data found here indicate that the RNA aptamer employs a conformational selection mechanism for binding theophylline in the absence of Mg(2+). The kinetic data found here also explain a very unusual property of this RNA-theophylline system: slow exchange on the NMR chemical shift time scale for a weakly binding complex. To our knowledge, it is unprecedented to have such a weakly binding complex (K(d) approximately 3.0 mM at 15 degrees C) show slow exchange on the NMR chemical shift time scale, but the results clearly demonstrate that slow exchange and weak binding are readily rationalized by a small k(on). Comparisons with other ligand-receptor interactions are presented.

The cavity-chaperone Skp protects its substrate from aggregation but allows independent folding of substrate domains.

Outer membrane proteins (OMPs) of gram-negative bacteria are synthesized in the cytosol and must cross the periplasm before insertion into the outer membrane. The 17-kDa protein (Skp) is a periplasmic chaperone that assists the folding and insertion of many OMPs, including OmpA, a model OMP with a membrane embedded beta-barrel domain and a periplasmic alphabeta domain. Structurally, Skp belongs to a family of cavity-containing chaperones that bind their substrates in the cavity, protecting them from aggregation. However, some substrates, such as OmpA, exceed the capacity of the chaperone cavity, posing a mechanistic challenge. Here, we provide direct NMR evidence that, while bound to Skp, the beta-barrel domain of OmpA is maintained in an unfolded state, whereas the periplasmic domain is folded in its native conformation. Complementary cross-linking and NMR relaxation experiments show that the OmpA beta-barrel is bound deep within the Skp cavity, whereas the folded periplasmic domain protrudes outside of the cavity where it tumbles independently from the rest of the complex. This domain-based chaperoning mechanism allows the transport of beta-barrels across the periplasm in an unfolded state, which may be important for efficient insertion into the outer membrane.

Measurement of imino 1H-1H residual dipolar couplings in RNA.

Imino (1)H-(15)N residual dipolar couplings (RDCs) provide additional structural information that complements standard (1)H-(1)H NOEs leading to improvements in both the local and global structure of RNAs. Here, we report measurement of imino (1)H-(1)H RDCs for the Iron Responsive Element (IRE) RNA and native E. coli tRNA(Val) using a BEST-Jcomp-HMQC2 experiment. (1)H-(1)H RDCs are observed between the imino protons in G-U wobble base pairs and between imino protons on neighboring base pairs in both RNAs. These imino (1)H-(1)H RDCs complement standard (1)H-(15)N RDCs because the (1)H-(1)H vectors generally point along the helical axis, roughly perpendicular to (1)H-(15)N RDCs. The use of longitudinal relaxation enhancement increased the signal-to-noise of the spectra by ~3.5-fold over the standard experiment. The ability to measure imino (1)H-(1)H RDCs offers a new restraint, which can be used in NMR domain orientation and structural studies of RNAs.

Solution structure of tRNAVal from refinement of homology model against residual dipolar coupling and SAXS data.

A procedure is presented for refinement of a homology model of E. coli tRNA(Val), originally based on the X-ray structure of yeast tRNA(Phe), using experimental residual dipolar coupling (RDC) and small angle X-ray scattering (SAXS) data. A spherical sampling algorithm is described for refinement against SAXS data that does not require a globbic approximation, which is particularly important for nucleic acids where such approximations are less appropriate. Substantially higher speed of the algorithm also makes its application favorable for proteins. In addition to the SAXS data, the structure refinement employed a sparse set of NMR data consisting of 24 imino N-H(N) RDCs measured with Pf1 phage alignment, and 20 imino N-H(N) RDCs obtained from magnetic field dependent alignment of tRNA(Val). The refinement strategy aims to largely retain the local geometry of the 58% identical tRNA(Phe) by ensuring that the atomic coordinates for short, overlapping segments of the ribose-phosphate backbone and the conserved base pairs remain close to those of the starting model. Local coordinate restraints are enforced using the non-crystallographic symmetry (NCS) term in the XPLOR-NIH or CNS software package, while still permitting modest movements of adjacent segments. The RDCs mainly drive the relative orientation of the helical arms, whereas the SAXS restraints ensure an overall molecular shape compatible with experimental scattering data. The resulting structure exhibits good cross-validation statistics (jack-knifed Q (free) = 14% for the Pf1 RDCs, compared to 25% for the starting model) and exhibits a larger angle between the two helical arms than observed in the X-ray structure of tRNA(Phe), in agreement with previous NMR-based tRNA(Val) models.

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme.

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca(2+), Mg(2+), Mn(2+), and Sr(2+) have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.

High-resolution pyrimidine- and ribose-specific 4D HCCH-COSY spectra of RNA using the filter diagonalization method.

The NMR spectra of nucleic acids suffer from severe peak overlap, which complicates resonance assignments. 4D NMR experiments can overcome much of the degeneracy in 2D and 3D spectra; however, the linear increase in acquisition time with each new dimension makes it impractical to acquire high-resolution 4D spectra using standard Fourier transform (FT) techniques. The filter diagonalization method (FDM) is a numerically efficient algorithm that fits the entire multi-dimensional time-domain data to a set of multi-dimensional oscillators. Selective 4D constant-time HCCH-COSY experiments that correlate the H5-C5-C6-H6 base spin systems of pyrimidines or the H1'-C1'-C2'-H2' spin systems of ribose sugars were acquired on the (13)C-labeled iron responsive element (IRE) RNA. FDM-processing of these 4D experiments recorded with only 8 complex points in the indirect dimensions showed superior spectral resolution than FT-processed spectra. Practical aspects of obtaining optimal FDM-processed spectra are discussed. The results here demonstrate that FDM-processing can be used to obtain high-resolution 4D spectra on a medium sized RNA in a fraction of the acquisition time normally required for high-resolution, high-dimensional spectra.

Imino proton exchange rates imply an induced-fit binding mechanism for the VEGF165-targeting aptamer, Macugen.

The 2'-fluoro/2'-O-methyl modified RNA aptamer Macugen is a potent inhibitor of the angiogenic regulatory protein, VEGF165. Macugen binds with high affinity to the heparin-binding domain (HBD) of VEGF165. Hydrogen exchange rates of the imino protons were measured for free Macugen and Macugen bound to the HBD or full-length VEGF to better understand the mechanism for high affinity binding. The results here show that the internal loop and hairpin loop of Macugen are highly dynamic in the free state and are greatly stabilized and/or protected from solvent upon protein binding.

Comparison of alignment tensors generated for native tRNA(Val) using magnetic fields and liquid crystalline media.

Residual dipolar couplings (RDCs) complement standard NOE distance and J-coupling torsion angle data to improve the local and global structure of biomolecules in solution. One powerful application of RDCs is for domain orientation studies, which are especially valuable for structural studies of nucleic acids, where the local structure of a double helix is readily modeled and the orientations of the helical domains can then be determined from RDC data. However, RDCs obtained from only one alignment media generally result in degenerate solutions for the orientation of multiple domains. In protein systems, different alignment media are typically used to eliminate this orientational degeneracy, where the combination of RDCs from two (or more) independent alignment tensors can be used to overcome this degeneracy. It is demonstrated here for native E. coli tRNA(Val) that many of the commonly used liquid crystalline alignment media result in very similar alignment tensors, which do not eliminate the 4-fold degeneracy for orienting the two helical domains in tRNA. The intrinsic magnetic susceptibility anisotropy (MSA) of the nucleobases in tRNA(Val) was also used to obtain RDCs for magnetic alignment at 800 and 900 MHz. While these RDCs yield a different alignment tensor, the specific orientation of this tensor combined with the high rhombicity for the tensors in the liquid crystalline media only eliminates two of the four degenerate orientations for tRNA(Val). Simulations are used to show that, in optimal cases, the combination of RDCs obtained from liquid crystalline medium and MSA-induced alignment can be used to obtain a unique orientation for the two helical domains in tRNA(Val).